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Abstract: Biorthogonal two-direction wavelet packets with dilation factor are introduced and their 
properties are discussed by means of the matrix theory and operator theory. A new ap- 
proach for constructing biorthogonal two-direction wavelet packets is developed. The for- 
mulae for performing iteration and decomposition are established. New Riesz bases for 
L? (R) are obtained by the given biorthogonal two-direction wavelet packets. Finally, an 
example for constructing biorthogonal two-direction wavelet packets is given. 
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1 Introduction 


Wavelet transform is an important mathematical tool with which data or functions can be 
divided into different frequency components. It is well known that there is a limitation for the 
time-frequency localization of a single wavelet. In other words, an orthogonal wavelet function 
with compact support and certain regularity is not symmetricl!:?]. Therefore, Geronimo et all’ 
constructed two functions ~(x) and wW2(x) whose translations and dilations form an orthonor- 
mal basis in L?(IR). The importance for these two functions is that they are continuous, well 
time-localized (or short support), and symmetry. This example tells us that if multiwavelets 
are used in an expansion, then better properties can be achieved. 

Orthogonal wavelet packets were firstly introduced by Coifman and Meyer'4). Orthogonal 
wavelet packets are used to further decompose wavelet components. Wavelet packets, due to 
their nice characteristics, have been widely applied to signal processing!*), image processing!®) 
and so on. Biorthogonal wavelet packets were given by Daubechies and Cohen!”). Biorthogonal 
wavelet packets are more flexible in application. In addition, wavelet packets provide better 
frequency localization than wavelets while time-domain localization is not lost. 

The two-scale refinable equation with scale a (2 < a € Z) 


p(t) = $ puplat — u) (1) 
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plays a basic role in the construction and application of scalar wavelets'®9l. Yangl!°! generalized 
the two-scale refinable equation, and established the biorthogonality criteria for two-direction 
refinable function and two-direction wavelets. Motivated by [10,11], we give the definition of 
biorthogonal two-direction wavelet packets and discuss their properties by means of matrix 
theory. The formulae for performing iterations and decomposition are also established. 


2 Two-direction multiresolution analysis with dilation a 


We begin with recalling some basic notations and results used later. Let R and C be the sets 
of all real and complex numbers, respectively. Denote by Z and Z+ the set of all integers and 
nonnegative integers, respectively. To obtain a uniform method for constructing biorthogonal 
two-direction wavelet with dilation a, let us give two-direction multiresolution analysis. 

Define the operators R, T and D as follows: 


(RAW = ft) (TAE) =ft-1), (Daft=ai flat), VfELR), tER. 


Then R, T and D are unitary operators on the Hilbert space L?(R). 
Definition 2.1 Given 2 < a € Z and ¢ € L?(R). If there exist {p} }ucz € (Z) and 
{pz Juez € €2(Z) such that 
=> pla DaT" + > pia? DaT "RG, (2) 
ucz uez 
then ¢ is said to be a two-direction refinable function (TDRF). The sequences {p} }ucz and 
{pz }uez are called the positive-direction mask (PDM) and negative-direction mask (NDM) of 
¢, respectively. 
The equation (2) can be simply written as 
olt) = J pi gat - u) + > pu plu- at); (3) 
ucz uEZ 
which becomes two-scale equation (1) in the case of p, = 0. 
By taking the Fourier transformation for the both sides of (2), we have 


(w) = pt (e) (w/a) +p (ce) (w/a), (4) 


where 
pt(z)=(1/a) $ ptz", z=e= 


ueZ 
is called the positive-direction mask symbol (PDMS) and 
= 1 = 
P (z) z a Yr z“ 
ucz 
is called the negative-direction mask symbol (NDMS). 


Proposition 2.1 Let ¢ be a TDRF. If there exist sequences {p} }uez and {py }uez of 4, 
then @ satisfies the following equation 





(w) pt(ea) p- (e™) | | d(w/a) (5) 
plu) ptle™) pole) || w/a) 
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where pt (z) is the PDMS and p7 (z) is the NDMS. 
Proof By taking the Fourier translation, we obtain (4). On the other hand, we rewrite 
(2) as 


Rọ =X phar? DTR +Y pean? DaT “gp. (6) 
ueZ ucz 
Simply 
= So pid(—at —u) + So puolu + at). (7) 
ueZ uceZ 
Also, by implementing the Fourier transformation for the both sides of (6), we have 
Jw) = pt (ee) pwa) + p-(e=) $ (w/a). (8) 


(4) and (8) can be written as (5). The proof is completed. 

Let ¢ be a two-direction refinable function with masks {p} }uez and {p} }uez. By virtue 
of the positive-direction mask {p{} ez and the negative-direction mask {pz }uez, we construct 
the following matrix equation 


¢(—t) PIu PË, 
a(t) = = B(at — u). (9) 
| o(t) | =| Pù Pu | 


It is easy to see that (5) and (9) are equivalent. The matrix 


1 Ia Pt = 
P(z)=->> amet Zz’, z=e7, 
eZ pt Pu 

is called the matrix mask symbol of ®(t). 

Definition 2.2 A pair {¢, $} of two-direction refinable functions is said to be biorthogonal 
if 

(¢, T*¢) = 50,k ($, T“R¢ ) =0, V k, uE Z, (10) 

where ĝo, is the Kronecker symbol. 

For a function ¢ € L?(R), we define a subspace sequence V; € L?(R) b 


V; = V ({DiT*9},e2 U {DAT" Ro} ez )- (11) 


Definition 2.3 The sequence {V;}jez defined by (11) is said to be a two-direction mul- 
tiresolution analysis (TDMRA) with scale a generated by ¢ if the following conditions are 
satisfied. 

(i) Vn C Vp4i, for all n € Z; 

(ii) AjezV; = {0}, Ujez V; is dense in L?(R); 

(iii) J € Vo if and only if Dif € V}, for all j € Z; 

(iv) The sequence {T*¢}xez U{T” Ro} nez is a Riesz basis of Vo. 

Let A = {1,2,---,a@— 1}. 

Definition 2.4 Let {V;}j;ez be a TDMRA generated by ¢. A subset {y, : 4 € A} of L?(R) 
is called a two-direction wavelet set (TDWS) associated to ¢ with scale a if 


{T" h. rez U {T" Ry, hez LEÀ, 
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forms a Riesz basis for Wọ where Wo = Vi © Vo. 

Proposition 2.2 Let {, :ı € A} of L?(R) be a TDWS associated to ¢ with scale a. If 
there exist sequences {q}, }uez and {q,,}uez of {wy : t € A}, then {y, : 4 € A} satisfies the 
following equation 





z Sy iw Fp TE 
T, (w) = oe a ih ee) aS a oa) , LEA, (12) 
p(w) qi(e=) gr (e™) (2) 


where 


u = 1 = Zii 
qt (2) = (1/a) J ai", @e==digz", zme. 


ueZ ueZ 
Proof For {y, :ı € A}, we have 


he => at o- P DeT*6+ Yo q,a DaT “RO, LEA, (13) 
ueZ ucz 
This equation can be simply written as 
w(t) =Y at plat-u) +) az, o(u-at), LEA. (14) 
ueZ ucz 


Implementing the Fourier transformation for the both sides of (13) yields 


bw) = qt (0 )dw/a) + ar (e™) w/a), LE. (15) 
We rewrite equation (14) as 
hht) =J plat- u) + Dia outat), cea. (16) 
ueZ ucZ 


The refinement equations (14) and.(16) lead to the following formula 


¥.(-t) qu, qtu t 
Y, (t) = = : " |@(at—u), LEÀ. (17) 
| p(t) | x | Gir ae 


' 


The frequency field form of the relation formula (17) is (12). The proof is completed. 


3 Two-direction wavelet packets 


We shall introduce the two-direction wavelet packets and investigate their properties. Let 
ġ and ¢ be a pair of biorthogonal two-direction refinable functions, we rewrite the symbols as 


go(t)=4(t), pL =p}, O =p, O= PO =o LEA, wed, 


dolt) =0t), PO =pt, PrO =p, PBO=Gt, POS LEA, ved, 


1 pa pa 1 pe Du A) 

(a) = — 7 a P(A) ae fe u = —iw/a 

P`“ (z) = ee ie koe | PX’ (z) D US. ao as , tea. 
ueZ | Pia Pu ueZ | Diu —u 
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Definition 3.1 Two families of two-direction refinable functions {Wan+a(t), n € Zt, à= 
0,1,---,a@— 1} and {hon+xr(t), n € Z+, à = 0,1,--- ,a— 1} are called biorthogonal two- 
direction wavelet packets (TDWPs) with respect to a pair of biorthogonal two-direction scaling 
functions @(t) and (t), respectively, if 


Pan+a (t) = Y pt dalat —u)+ X pr Yn (u = at), (18) 
ueZ uez 

Daner(t) = Y pt dn(at — u) + X pu alu — at). (19) 
ued ueZ 


By implementing the Fourier transformation for the both side of (18) and (19), respectively, we 


have 
Dans (w) = P by (w/a) + PO Yn (w/a), A=0,1,---,4-1, (20) 
fanla) = Pt fn (w/a) + P-M dn w/a), A=O A -1, (21) 
where 
1 f 1 i 
PHO (wy) = = HA gine -Awu = — —(A) grin er deere, 
(w) DD emne Pw) A iu et’ > i=0,1,---,a—1. (22) 
ueZ ucz 


~ 1 ; Z 1 ; 
PHN) (w) E pe oe em, P- w) Dea Ao p enri Le Odie a 1. (23) 
a a 
ucz uez 
Properties and advantages of the biorthogonal two-direction wavelet packets with dilation a are 
investigated as follows. By applying the same method, one can obtain orthogonal two-direction 
packets with dilation a given by (18). 
Lemma 3.1 Suppose ®(t) and &(t) are biorthogonal scaling function vectors. P(z) and 


P(z) are their matrix symbols, respectively, and wj (j = 1,2,---,@) are roots of equation 
z” —1=0. Then 
Y Phwjz) Pwj)" =, 2-1 = 0. (24) 
j=l 
It is equivalent to 
Y PiPitak = aôo,klr, k €Z. (25) 
LEZ 


We can easily prove Theorem 3.1 by applying the lemma. 

Theorem 3.1 Suppose that y, (t) and ,(t) are biorthogonal two-direction wavelets asso- 
ciated to $(t) and ġ(t), respectively, wj (j = 1,2,- +> ,a) are roots of equations z* — 1 = 0 with 
z=e~*. Then 


a 
Y PO (wz) Pwj) =O, LEA, (26) 
j=1 
Y P(wjz) PO (wz) =O, LEA, (27) 
j=l 
Z PO (wz) PO (w;z)* = a-i LEA (28) 


j=1 
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Theorem 3.2 Suppose n € Z and 


n=)  Bjol*, B; € {0,1,2,--- a — 1}. (29) 
j=l 
Then 
dn(w) = [| PO) (e77 )bo(0), (30) 
j=l 
Dnlw) = J| PO (c= )o(0). (31) 
j=1 


Proof When n = 0, (30) follows from (4). Suppose (30) holds for 0 < n < a”, When 
a’ < n < atl, by using (20) and the inductive assumption, we have 


Ûn(w) = PO) (c=) bray 


co lee} 


PPD) (e) II PP) (ez7FT Jbo(0) = Il Pi) (eS )yho(0). 


j=l j=l 


Then (30) holds by induction. We can also obtain (31) by using the same argument. The proof 
is completed. 
Theorem 3.3 Suppose k, l € Z and n € Z,. Then 


(bn(z Tz k), Pn(z e; 1) r ôk Ir. (32) 


Proof When n = 0, (32) holds. Suppose (32) holds for 0 < n < a7. When a < n < 
a’+1, by using the inductive assumption (26), we have 


(dn (t mE k), Vnl(t m l)) 


= È f iontas 


1 2an(j+1) p 2 i 
= J Po) (2) POD (2) bragged 
jez 2anj g 


2am ie i 
ee PO)(z) POI*(z) Y > hay(w + 2amj) - Vpaj(w + 2arjje™ -edw 
on 0 jez Q a 
207 
= L [7 pz) PO (z)e-ik-De guy 
2r Jo 
1 2r Od 
= Ss (81) r>(G1)* —i(k—-lw q 
Zn f ( oF (Wmz)P (wmz) Je w 
1 


2r , 
= — Pain i oe dw = ôk Ir. 
2a 0 i 
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Then formula (32) is established by the inductive assumption. The proof is completed. 
Theorem 3.4 Suppose k, l € Zand n € Z and \ € {1,2,--- ,a—1}. We have 


(han( — k), Panta(t — l)) = 0. (33) 
Proof By applying (18), (23), we have 
(hanla E k), Panya lT = k)) 
= 5 f POOP inuia (ade Mau 


2an i l 
= x P (2) PO*(z) > Yn(w + 2mar) + Yn (w + 2mar)je™ -Ddw 
T Jo 
mEZ 
1 2an r , 
= a PA (z) P91*(z)eHA—De au) 
T Jo 
1 2an a-—1 N l 
= ( Pp) (wy z) PO* (wz) Je ied =0. 


Qn 
0 P= 


Therefore, for any k, l € Z, (33) is established. The proof is completed. 
4 The direct decomposition for space L?(R) 


In this section, we will decompose subspace Vj, V; and W;, W, by virtue of a series of two- 
direction wavelet packet subspaces. Furthermore, we present the direct decomposition for space 
L?(R). Define 


Un = V ({1*bn}rez (PR when): 
Ün = V (11 Pa}rez U {TRE hez) 


and define a dilation operator (S¢)(t) = ¢(at) where (t) € L?(R). We have 
S'Un = Vain ® Ugins1 ®Uainzat-1, l nE Z+, 
S'Û, = Tain Ð Going ® Ü inkai l, ne zt 


and S°=1, S! = S,--- , SÌL. The following result is equivalent to Theorem 3.4. 
Lemma 4.1 Let {yanza(t), n € Zt, A = 0,1,--- ,a— 1} and {¥anza(t),n € Zt, A = 
0,1,-+- ,@—1} be biorthogonal TDWPs. Then we have the following decomposition formulus: 


a~l 


1 ry (7}* 
plat — k) = aa > PBD banil —1), kez, (34) 
j=1 1EZ 
? HE PORE 
plat- k) = -39 9 Putant- l), kez. (35) 


j=1 leZ 
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Proof By the definition of TDWPs and applying (24), we have 


a-i 
ZED anl- i) 


j=1 leZ 
Law 2) 
= POLY PP yalt- ad — m) 
jJ=1 led meZ 
i Be ha 
a2 3 2 2 POr PO app (crt —al—m) 
j=1 IEZ mEZ 
KE PS 
æ S D D POPO bs a(at — 8) 
j=1 lEZ sEZ 
ie ait . 
= a yD heo y Fa Pe, 
j=1 sEZ leZ 
= Palat — k). 


Similarly, one can obtain (35) by using the same method. The proof is completed. 


VOL. 27 


Lemma 4.2 Let {Yan+ar(t), ne Z+, à= 0, 1,--- a — 1} and {Vana (t), ne Zt, A= 


0,1,--- ,@—1} be biorthogonal TDWPs. Then 


(Um(t ~ k), Yn(t — l)) = bm nk, Ly, m, n= 0,1,--- : 


Theorem 4.1 Let j =0,1,---, and 

Vi = ®osn<aiUn, Wy = ®o<ncaiUn, 
Z = @o<n<aiUn, wW; E Po<n<os Un, 
where “@®” denotes direct sum of subspaces. Then 


L? (R) = @rnezUn = @nezUn, 


where 
UmLUn, m,n€Zy, mén, 


k, LEZ. 


(36) 


(37) 


(38) 


(39) 


(40) 


and {m,(t—k), Ym, (k — t), i = 1,2, k,l € Z}, {tn,(t —k), Ùn (k- t), = 1,2, k,l € Z} are 


biorthogonal bases of U, and U,, respectively, where m, n € Z4. 


Proof Noting that Vo = Up, W = @%,U, and V; = SIWo, we can obtain (37). Since 


{V;};ez and {V;},;ez are all MRAs of L?(R), one obtains 


L?(R) = Vo @ ( @jez Wj) = Vo @ ( Byez W,). 


This completes the proof of Theorem 4.1. 


NO.5 Li Lan et al: Biorthogonal Two-direction Wavelet Packets with a Positive Integer Dilation Factor 909 





5 An Example 


In this section, we give an example to demonstrate the general theory in section 3. 
Let (t), o(t) be a pair of biorthogonal two direction refinable functions with dilation factor 


2. If they satisfy 
Olt) = SAt + 1) + (24) + ZARE- 1) + A2 +1) - Eol = 1), 
a(t) = Sga +1) + (2t) + igot —1)+ “Z-z yije “Z-z -1), 
and 
o(-t) = —29(—2¢ + 1) + Z624) — Zo(-2¢ — 1) — At + 1) + 42t — 1), 
Bt) = -Spt 1) + EB) + EG-t- 1) - Et 1) + Goce), 


then y(t), w(t) are biorthogonal two direction wavelets associate to ¢(t), o(t), respectively, 
which satisfy the following equations 


p(t) = poet +1) — o(2t) + Solat —1)+ zoz +1)— Toz -—1), 


2 Igor +1) —G(2t) + Tg api EZ-a ie 53-2 cii 





and 
sej=- otar- Lota Aoa y - Moar + 1) o-i, 
o( fs on me o( an +1) — ig an) — SVT §( #-1)- Garr 1+ 2 gtt- 2): 


By (30)-(35), we can obtain the biorthogonal two-direction wavelet packets. 
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